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Abstract 

The vibrational properties of a high-quality realistic model of amorphous silicon 
are examined. The longitudinal and transverse dynamical structure factors are cal- 
culated, and fitted to a damped harmonic oscillator (DHO) function. The width 
r of the best-fit DHO to the longitudinal dynamical structure factor scales ap- 
proximately as k 2 for wavevectors k < 0.55A , which is above the loffe-Regel 
crossover frequency separating the propagating and diffusing regimes, occurring at 
k = 0.38±0.03A . Using the DHO function fitting function for the transverse 
dynamical structure factor (without theoretical justification), gives a dependence of 
r oc k a with a ~ 2.5 for wavevectors k < 0.7A . There was no evidence for r oc k 4 
behaviour for either polarization. 
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1 Introduction 



The physics of the atomic vibrations present in amorphous materials has been 
the subject of much investigation, see e.g. ([H; Q; 0). Advances in experimen- 
tal techniques such as inelastic neutron scattering (INS) (0), inelastic X-ray 
scattering (IXS) (0), hyper- Raman scattering (0) and terahertz time-domain 
absorption spectroscopy (jf!) have allowed vibrations with frequencies of the 
order of one terahertz (1 THz ~ 4.1 meV/h) to be examined. Vibrations at 
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these frequencies are of particular interest, as these frequencies are compara- 
ble to those of the boson peak (an excess of vibrational modes as compared 
to the Debye elastic waves) M) and of the Ioffe-Regel (IR) crossover between 
propagation and diffusion (0; @; @) observed in amorphous materials. 

Vibrations in amorphous materials can be characterised by the spectral den- 
sity. The spectral-density operator is A(e) = S(e — D), where e = uj 2 is 



the squared vibrational frequency, and D is the dynamical matrix (llOl ). In 
the plane-wave basis, the matrix elements of the spectral-density operator 
A^(e) = (k/3|A|k/3), where k is the wavevector and /3 is the polarization of 
the plane wave, are proportional to the dynamical structure factor Sp(k, uS) 
(0), which, for the longitudinal polarization (/3 = L), is the quantity measured 
by INS or IXS (0). 

This proportionality is straightforward to see. Into the expression for Aj C| g(e) 
above, the resolution of the identity in the basis of the eigenvectors { |e^) } of 
the dynamical matrix D is substituted to give: 

A k/3 (e) = E(k/3|e l )(e J |5( e -E))|e J )(e J |k/5), (1) 



where i and j are labels for the eigenvectors, which are normalized to unity, 
(ej|ej) = 8ij. This expression is equal to 

A^eHEKe.lk/^e-e,), (2) 



where q is the energy of eigenmode i. Transforming equation [2] into frequency 
space, via e = u 2 , gives 

A ^M = ^EK^I k ^)l 2 ^-^), (3) 



where u>i is the vibrational frequency of eigenmode i. Equation [3] is valid for 
a plane wave of any polarization (3. At high temperature, and in the single- 
excitation approximation, the dynamical structure factor for the longitudinal 



polarization, (3 — L, is defined as fllll): 



kmTk 2 

S 0=L (k,u;) = £ IfolW = L)\ 2 5(uj - (4) 



where ks is Boltzmann's constant, T is the temperature, and m is the atomic 
mass, and the Debye- Waller factor has been ignored (as in (7)). This is clearly 
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proportional to the spectral-density operator A^(uj) in equation [3j with po- 
larization {3 = L. By defining a plane wave as 



M> = 7^E n k /^ ik - rj , (5) 



where N is the number of atoms, is the polarization vector of the plane 
wave, and Tj is the position vector of atom j, and substituting this into equa- 
tion HI an expression for the dynamical structure factor for any polarization 
j3 can be derived as 

S,(k,u) = k Jn? r Y. I E( n ^ • e? } )e- ik - r f 8{u> - m), (6) 
mar ~r , 



where is the displacement vector of atom j vibrating in eigenmode i. The 
transverse dynamical structure factor 5r(k, uS) is not readily accessible in scat- 
tering experiments, but can easily be calculated from simulation. In practice, 
the orientationally averaged dynamical structure factor Si(/c, o>) is obtained 
from scattering experiments (|3j), and, in this paper, both the orientationally 
averaged dynamical structure factors Si t T(k,u>) are calculated for a-Si. 

As observed in experiment, Si,(k,u) of an amorphous material exhibits a sin- 
gle peak for low values of wavevector k (0; S)- The frequency at which this 
peak is located increases linearly with increasing k (at least for low k); this 
is the regime of linear dispersion. The width of this peak also increases with 
increasing k, but more rapidly. Eventually, the width of the peak therefore 
becomes comparable to the peak frequency. The frequency at which this hap- 
pens (typically ~ 1 THz) is called the Ioffe-Regel (IR) crossover frequency 



(1121 ). Above this frequency, the mean free path of the wave is shorter than 



the wavelength, and the external waves do not propagate in the conventional 



sense, instead transferring energy by diffusion (jl3l ). 



These basic features are undisputed, but large uncertainties in the experimen- 
tal data (see, e.g., (0)) have meant that it has been difficult to decide between 
the two main theoretical models for the behaviour of vibrational plane waves 
in amorphous materials. The first of these (0) suggests that the memory func- 
tion is the sum of two processes: structural relaxation, assumed to be too slow 
to have an effect on vibrational timescales, and an instantaneous (or at least 
very fast) process, usually represented as a <5-function in time (0). This leads to 
Sz,(k, co) having the form of a damped harmonic oscillator (DHO) function (Q). 
Fitting a DHO function to the experimentally observed peak gives a width T 
which scales roughly as the square of the wavevector, r oc k 2 , for many mate- 
rials p). This relationship is often seen to hold up to frequencies above the IR 
crossover. The second of these models (see, e.g., (Il5l)) is based on an effective- 
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medium approximation (EMA) (Il6l ) and asserts that, in the region of the IR 
crossover, a DHO function is not a valid fit to Si,(k,u}), and that the scatter- 
ing becomes much stronger (see, e.g., (fla)). A functional form for S^kjUj) is 
used which imposes r oc k A by construction, and this has given good fits of 
the dynamical structure factor at frequencies close to the IR crossover (114l ). 
The behaviour r oc k A was also found recently on fitting a DHO function to 
Sl(/c, u>) of lithium diborate glass (|2j). Recent inelastic ultra-violet scattering 



measurements (1171 ) on vitreous silica suggest that there may be three regimes 



at different frequencies, with r oc k 2 , T oc k 4 and r oc k 2 respectively. 

The simplest possible model for the effect of disorder on vibrations, a force- 
constant-disordered lattice, gives T oc k 4 (jUjEi). A theoretical model 00) 



which includes positional disorder predicts r oc k 2 , but the model is one in 
which the dynamical matrix is assumed to be a Euclidean random matrix, and 
the positions of the atoms are taken to be uncorrelated and random. It is not 
clear how representative such a model is of real amorphous systems. 

Further investigation of the behaviour of the width T(k) of Si{k,uj) might 
therefore be useful to decide on the generality or otherwise of these models. 
Investigating S^k^uj) by simulation is helpful, because the large error bars 
which characterize the experimental results are absent. Simulated vitreous sil- 
ica has been shown to have r oc k 2 (jsl; 22) using a DHO function, and fitting 



Lorentzian functions to Si(k, u) of amorphous silicon (a-Si), simulated with 
the Stillinger- Weber potential (1231 ). also resulted in a T oc k 2 dependence for 
k < 0.3A~ (24). In this paper, we also examine the wavevector dependence 
of T in a simulated model of a-Si, but one simulated using the more phys- 
ically realistic modified Stillinger- Weber potential (jiil ). as described in the 
next section, as well as fitting to S^k, uj) using the DHO fitting function. Fit- 
ting Lorentzian functions to S^k^uj) of a-Si, produced behaviour suggestive 
of a T oc k 2 dependence for k < 0.3A (1261 ) . but this model was also con- 
structed with the Stillinger- Weber potential. The behaviour of both Si,(k,uj) 
and ST{k,u) is investigated here up to wavevectors of k ~ 2.0A , which 
allows us to show where the DHO fitting function breaks down, and to inves- 
tigate the vibrational properties of this model of a-Si at frequencies above the 
Ioffe-Regel crossover. 



2 Methods 



The model of a-Si studied here is a 4096-atom cubic model with box-size 
43.274A created by G. T. Barkema using a modified WWW method (j27|), 
and then relaxed by us (1281 ) using the modified Stillinger- Weber (mSW) inter- 
atomic potential (1251 ). The dynamical matrix was diagonalized using a Lanczos 
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routine (1291 ). The dynamical structure factors were calculated using equation [6j 
To remove unimportant physical constants, the modified dynamical structure 
factor S%(k,u>) = -^,Sp{k, uj) is plotted here, where the polarization (3 = L or 
T, as in These modified dynamical structure factors were fitted either to 
a single DHO function with two free parameters: peak frequency Q and width 
T, or to a sum of two DHO functions (for certain wavevectors) with five free 
parameters: peak frequencies and widths of the two DHO functions, and the 
ratio of their intensities. The peak normalization is fixed via / Ak^(e)(ie = 1. 
To prevent artefacts, all data below uj = 2.5 meV were ignored in the fitting; 
the lowest non-zero vibrational frequency for this model is uj ~ 3.7 meV (at 
this frequency, the modes are of mostly transverse character, as seen by the 
peak in Srik^uj) at this frequency in Figured]). Due to the finite size of the 
model, the lowest value of wavevector k compatible with the periodic boundary 
conditions is k = 0.145A . 



Strictly speaking, the dynamical structure factor is a weighted sum of delta- 
functions, with non-zero values only at frequencies equal to the mode frequen- 
cies (see equation E]) • However, for ease of representation, the delta- functions 
in the spectral densities were broadened, before fitting, by a Lorentzian func- 
tion with width r br . The width r br of the broadening chosen has a small, but 
noticeable, effect on the value of the width T of the best-fit DHO function, 
and hence on the value of the exponent a in the expression r oc k a . r br was 
usually chosen to be larger than a few energy-level spacings (the mean energy- 
level spacing in this model is ~ 0.007 meV), to smooth the data effectively for 
fitting. However, the larger r br becomes, the greater is the effect on the raw 
data. 



Increasing the value of r br tends to decrease the value of a. This is to be 
expected, if we assume that the effect of the Lorentzian broadening is to 
increase the true width by an amount comparable to r br . If the true widths 
are plotted on a double-logarithmic plot, then adding a constant T br to each 
width will shift the low-r points further up the plot than the high-r points. 
This will have the overall effect of reducing the gradient of the straight-line 
fitted to the broadened data. A value of r br is required which provides a 
compromise between producing broadened data that are smooth enough to be 
fitted reliably, but giving small enough broadening so that the data are not 
greatly affected. For the longitudinal polarization, values of r br in the range 
0.03 — 0.07 meV were used (the narrowest peak has width ~ 0.2 meV), and 
for the transverse polarization, values of r br in the range 0.015 — 0.04 meV 
were used (the narrowest peak has width ~ 0.15 meV). 
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Fig. 1. The modified longitudinal dynamical structure factor S^(k,u) of the realistic 
a-Si model, for values of wavevector k < 0.6A , and the corresponding best-fit 
DHO functions. The width Tb r of the broadening used was 0.05 meV. 

3 Results 



Longitudinal polarization 

Figure [I] shows representative longitudinal modified dynamical structure fac- 
tors Sl(k, u) for values of wavevector k < 0.6A , together with their best-fit 
DHO functions. It is clear that, in this wavevector range, S* L (k,u) has only 
one peak, and that, for these wavevectors, the dynamical structure factor is 
well described by the DHO function. The DHO function sometimes deviates 
slightly from the data, e.g. on the low-energy side at around to ~ 10 meV of 
the k = 0.290A peak, but these deviations are small. 

The width Y of the best-fit DHO function is presented in Figure H as a 
function of wavevector k, on a double-logarithmic scale, for a representa- 
tive value of r br = 0.05 meV. In Figure [2J the data points are seen to lie 
roughly on a straight line, indicating that a power-law dependence, Y oc k a , 
is an appropriate description. Other values of Y^ T give a width Y which ex- 
hibits power-law behaviour as well. Figure [2] also shows the best-fit Y oc k 2 
(a = 2), T oc k A (a = 4), and unconstrained a lines to the data in the range 
0.145 < k < 0.5A~ . If a is allowed to vary, then the best-fit value of a for 
these data is 2.24 ± 0.07 (this error is estimated by using the sum-of-squares 
deviation of the data from the best-fit straight line). The best-fit Y oc k 2 line 
is shown to be close to the data, while Y oc k 4 does not represent the data at 
all. 
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Fig. 2. The width T of the best-fit DHO function to the modified longitudinal 
dynamical structure factor S^(k, u>) of the realistic a-Si model with Tb r = 0.05 meV, 
on a double-logarithmic scale, with best-fit a = 2, a = 4, and a-unconstrained lines, 
where r oc k a , fitted to the data in the range 0.145 < k < 0.5A . 
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The precise value of a obtained from data such as those in Figure [2] depends on 
the value of broadening r br used, and the range of wavevectors studied. It will 
become clear from the next section (Figure H]) that the power-law behaviour 
r oc k a is not obeyed for k > 0.55A . Hence data sets were generated for 
values of broadening 0.03 < Tb r < 0.07 meV, and ranges of wavevector start- 
ing at 0.145 < k min < 0.26A (to remove possible finite-size effects), and 
extending up to 0.45 < k max < 0.55A . For these data sets, the mean value 
of a was 2.08, with a standard deviation of 0.15. 98% of the data sets had a 
within ±0.3 of this mean. For no value of broadening Tb r or wavevector range 
was a line with a = 4 found to be close to the data. 



The Ioffe-Regel (IR) crossover for this model of a-Si was estimated, defining 
it to occur at T ~ Vl/2-k (J3), where Vt is the excitation frequency in S^(k,cu). 
For a given value of and range of wavevector k, an estimate of the error in 
the IR crossover frequency was made, by using the sum-of-squares deviation 
from the best-fit straight-line (on a double-logarithmic plot) to obtain errors 
on the value of a, and the sound velocity. As a first approximation, the IR 
crossover was assumed to lie in the centre of the overlap region of the two 
sets of error bars. Figure [3] displays this procedure for r br = 0.05 meV and 
wavevectors 0.145 < k < 0.5A . The longitudinal IR crossover frequency 
was found to be Qirx = 19.1 ± 1.3 meV, corresponding to a wavevector of 
km,L — 0.38 ± 0.03A . This is somewhat lower than the previous estimate 



of ~ 25 meV (124 ). due, presumably, to the different interatomic potential and 
fitting function used here. Different values of Tb r and wavevector range give 
values of SIir } l within the error bars given above. 

It follows from this value of £Iir,l and Figure [3] that the approximate r oc k 2 
behaviour for longitudinal modes for this model of a-Si continues to frequencies 
and wavevectors above the IR crossover frequency and wavevector, as found in 
experiment and simulations for other materials (J3j), if the DHO fitting function 
is used. There is no sharp change in the behaviour of S* L (k,u) evident at or 
around the IR crossover frequency for this model (Figured]), nor in the pseudo- 
dispersion law fl = czk, with Cl ~ 7670 m s -1 (Figure [3]). 

At higher wavevectors, the DHO fitting function ceases to produce a good 
fit (Figure [5]). In addition, for wavevectors in the range 0.55 < k < 0.8 A , 
the peak in S1(k,u) narrows (Figure Hj), although the excitation frequency 
continues to increase monotonically (Figure [5]). This narrowing is thought to 
be related to remnants of vibrational-mode localization occurring in the gap 
in the vibrational density of states between the acoustic and optic bands, at 
uj ~ 35 meV. (Peak excitation frequencies for 0.55 < k < 0.8A occur in 
this region.) This is shown in Figure [6], where the inverse participation ratio 
of an eigenmode, defined as (J2i e 2 ) 2 /NJ2i e t, where e« is the amplitude of the 
displacement of atom % in that mode, and N is the number of atoms, is shown 
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Fig. 3. The width T of the best-fit DHO to the modified longitudinal dynamical 
structure factor S* L (k : oj) of the realistic a-Si model, with broadening Tb r = 0.05 
meV, and the scaled peak frequency S1/27T, plotted against wavevector k, with errors 
to the best-fit lines. The IR crossover occurs in the region where the error bars 
overlap. 
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Fig. 4. The width T of the best-fit DHO to the modified longitudinal dynamical 
structure factor S* L {k,uj) of the realistic a-Si model with Tb r = 0.05 meV. 
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Fig. 5. The modified longitudinal dynamical structure factor S* L (k,uj) of the realistic 
a-Si model, for values of wavevector 0.6 < k < l.oA 1 , and the best-fit DHO 
functions. The width of the broadening used was Tb r = 0.05 meV. 




Fig. 6. The inverse participation ratio and vibrational density of states g(oj) of the 
realistic a-Si model. 
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and is seen to decrease (indicating greater localization) at roughly the same 
frequency (u ~ 35 meV) as the dip in the vibrational density of states between 
the overlapping bands. This anomalous dip in T(k) (Figure H]) at k ~ 0.7A is 
not associated with the IR crossover, which occurs at k IRL ~ 0.38A (Figure 
E]). 

At even higher wavevectors, k > 1.0 A , a second peak is evident in S^(k,uj) 
at low frequency (Figure ED- This peak is not accounted for either in the 
memory-function or EMA approaches, but in other materials is thought to 
be a signature of modes with a largely transverse character mixing with the 
longitudinal spectrum, this disorder-induced mode mixing being a result of the 
absence of purely polarized modes in an amorphous material, see, e.g., (J3 ; 19 ) . 



It will be clear from Figure [TTJ in the next section that the dominant peak in 
S^(k,uj) at high k indeed is responsible for the additional low-frequency peak 
in S1(k,uj). For sufficiently large values of wavevector k, the spectral density 
becomes similar to the VDOS (|7|). 

Figure [7] also shows the best-fit sum of two DHO functions, which give rea- 
sonable fits to the data. No quantitative analysis of the widths was performed 
for these wavevectors, as the memory-function approach does not predict the 
existence of this second peak, and so fitting the data to two DHO functions is 
largely done ad hoc. Also, for wavevectors 1.0 < k < 1.3A the second peak 
is small and flat. The widths found for this peak at these wavevectors vary 
significantly even for data from the same or very similar wavevectors, because 
the second-peak intensity is so small that its width can vary substantially 
without greatly affecting the sum-of-squares deviation. It is thus difficult to 
extract meaningful data about the peak widths for these wavevectors. 



Transverse polarization 



Figure [8] shows representative transverse modified dynamical structure factors 
S^(k,uo) J and their best-fit DHO functions, for wavevectors k < l.oA . The 
width Tbr of the broadening used was 0.02 meV. There is only one peak in 
S^(k,uj) at these frequencies. The memory-function approach does not pre- 
dict a particular functional form for the peak shape in the transverse spectral 
density, nor in fact does any current theory (30). Therefore, different func- 
tional forms for the fitting function were tried in an ad hoc fashion: Gaussian, 
Lorentzian and DHO. At very low wavevectors, k < 0.4A , the DHO and 
Lorentzian functions both gave good fits to the data. Above this wavevector, 
the DHO gave quantitatively the better fit to the data (in terms of having a 
lower sum-of-squares deviation), and this remained so for all wavevectors for 
which the data could reasonably be described by a single peak (k < 1.0A ). 
For this reason, and to facilitate a comparison with the data for the longitudi- 
nal polarization, only the results from the DHO fitting function are considered 
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CD (meV) 

Fig. 7. The modified longitudinal dynamical structure factor S* L (k,uj) of the realistic 
a-Si model, for values of wavevector k > 1.0 A , and the best-fit sum of two DHO 
functions. The value of broadening used was Tb r = 0.2 meV. 
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Fig. 8. The modified transverse dynamical structure factor S^(k,to) of the realistic 
a-Si model, for values of wavevector k < l.oA \ with the corresponding best-fit 
single DHO functions. The value of broadening used was Tb r = 0.02 meV. 

here. From Figure [U it is clear that the DHO function provides only a reason- 
able approximation to the data at low wavevectors. The noise in the data is 
due in part to the sparseness of eigenvalues at low frequencies. 

As for the longitudinal case, the precise value of the exponent a in the de- 
pendence r oc k a depends on the value of broadening r br used, and the range 
of wavevectors studied. For the transverse polarization, the physical relevance 
of the DHO best-fit width T is less clear, as there is no theoretical reason to 
expect the DHO function to be a good fit to S^(k,u). A discussion of such a 
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Fig. 9. The width T of the best-fit DHO to the modified transverse dynamical 
structure factor S^(k,tu) of the realistic a-Si model, on a double-logarithmic scale, 
with best-fit a = 2, a = 4 and unconstrained-a straight lines fitted to the range 
0.145 < k < 0.6 A . The value of broadening used was rv = 0.02 meV. 

fit is included here for completeness. Figure [9] shows the width T of the best- 
fit DHO functions to S^(k,u) as a function of wavevector k, as well as the 
best-fit lines with a = 2, a = 4 and unconstrained a, fitted to the wavevector 
range 0.145 < k < 0.6A -1 , and with T br = 0.02 meV. Above k ~ 0.6A~\ 
the DHO function gives a less good fit (Figure [8]), and deviations from linear 
dispersion are also present. The value of a in this case is 2.44 ± 0.05. It is 
clear that neither a = 2 or a = 4 represents the data well, although a = 2 
is closer. The precise value of a is likely to be strongly dependent on the fit- 
ting function used. For values of broadening 0.015 < Tb r < 0.04 meV and for 
ranges of wavevector starting at 0.145 < k m i n < 0.26A , and extending to 
0.5 < k max < 0.7A , the mean value of a found was 2.54, with a standard 
deviation of 0.14. 96% of the data sets generated had a within 0.3 of this 
mean. For no value of broadening or wavevector range did the exponent a = 4 
give a good representation of the data. 

Using the best-fit line for 0.145 < k < 0.6 A , the transverse IR crossover was 
estimated to be at a frequency of flm t r = 11.1 ± 0.5 meV, corresponding to a 
wavevector of kmp = 0.44±0.02A (Figure fTUj) . using the same method as for 
the longitudinal polarization. This is lower than the previous estimate of ~ 17 



meV (1241 ) . but the different interatomic potential and fitting function used 
in this study may explain this. As for the longitudinal polarization, no sharp 
change is evident in the width or the shape of S^(k, uS) at the transverse IR 
crossover (Figured]), nor in the pseudo-dispersion law fl = crk, with ct ~ 3820 
m s _1 (Figure fTUl) . 
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Fig. 10. The width T of the best-fit DHO to the modified transverse dynamical 
structure factor Sj,{k,w) of the realistic a-Si model, with broadening Tb r = 0.02 
meV, and the scaled peak frequency f2/27r, plotted against wavevector k, with errors 
to the best-fit lines. The IR crossover occurs in the region where the error bars 
overlap. 




Fig. 11. The modified transverse dynamical structure factor S^(k, ui) of the realistic 
a-Si model, for values of wavevector k > l.oA 1 . The value of broadening used was 
r br = 0.2 meV. 
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For wavevectors above k ~ 1.0 A , S^(k,uj) has one dominant peak at low 
frequency, with a flat shoulder at higher frequencies (Figure [TTj) . It is clear 
that a single peak fit is no longer an adequate description of the transverse 
dynamical structure factor in this wavevector region. Neither single DHO nor 
double DHO fits are very representative of the data for these wavevectors. As 
mentioned previously, the dominant low-frequency peak in S^(k, oS) appears in 
the longitudinal spectrum at high k (Figures [71TTT]) . and, vice versa, the high- 
frequency peak in S* L (k,u) also begins to appear in the transverse spectra, 
indicative of strong mode mixing at such high values of k. 



4 Discussion 



The width Y of the dynamical structure factor of amorphous materials in- 
creases with increasing wavevector k, but there is much dispute about the 
functional form of this increase. Broadly speaking, using a DHO fitting func- 
tion (as proposed in (3)) to fit to the dynamical structure factor peaks gives 
a behaviour of T oc k 2 at low k, with the exception of lithium diborate, for 
which T oc fc 4 (0) is found. A second method, based on an EMA approach 
(fl~6l). which involves modelling the dynamical structure factor with a function 
which has r oc k A by construction, has also yielded good fits at frequencies 
close to the Ioffe-Regel crossover between propagating and diffusive modes. 

This paper describes a model of amorphous silicon which is shown to exhibit 
an approximate r oc k 2 dependence for the longitudinal dynamical structure 
factor. This relationship holds for peak frequencies above the IR crossover 
frequency, as has been found for other materials modelled with the DHO 
function (]3j). An ad hoc fitting of the transverse dynamical structure factor 
to a DHO function (without theoretical justification) gives a slightly higher 
exponent, about 2.5. For neither polarization was any evidence for r oc k 4 
found, nor was there any sharp change in the behaviour of the dynamical 
structure factor at the IR crossover, as proposed in the EMA method. 

The nature of vibrations in amorphous systems in this frequency region is not 
yet fully understood. It is known that a force-constant-disordered crystalline 



lattice gives r oc k (Il8l ; Il9l ). while many positionally disordered materials, 



including a-Si, give F oc k , as does a positionally disordered analytical model 



(1201 ; |21| ). It may therefore be positional disorder which is responsible for the 
reduction in the exponent from 4 to 2, although the mechanism of this reduc- 
tion remains unclear. Of particular interest is the recent experimental work 
of Masciovecchio et al. (17), which suggests that there may be three regimes 
in vitreous silica, with respectively r oc k 2 , T oc k A and r oc k 2 . If this be- 
haviour is generally true, then these regions could occur at different ranges of 
wavevector for different materials, and this could account for different depen- 
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dences being observed between different materials, e.g. the r oc k 4 dependence 
in lithium diborate (|2j), and for the intermediate value of exponent a = 2.54 
found for the transverse polarization above. 
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